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Abstract —It is well known that spatially coupled (SC) codes 
with erasure-BP decoding have powerful error correcting capa¬ 
bility over memoryless erasure channels. However, the decoding 
performance of SC-codes significantly degrades when they are 
used over burst erasure channels. In this paper, we propose band 
splitting permutations (BSP) suitable for ( l,r,L ) SC-codes. The 
BSP splits a diagonal band in a base matrix into multiple bands 
in order to enhance the span of the stopping sets in the base 
matrix. As theoretical performance guarantees, lower and upper 
bounds on the maximal burst correctable length of the permuted 
(l, r, L ) SC-codes are presented. Those bounds indicate that the 
maximal correctable burst ratio of the permuted SC-codes is 
given by A max — I /k where k = r/l. This implies the asymptotic 
optimality of the permuted SC-codes in terms of burst erasure 
correction. 

I. Introduction 

Low-Density Parity-Check (LDPC) codes that are linear 
codes defined by extremely sparse parity check matrices were 
developed by Gallager in 1963 ||T). The combination of LDPC 
codes and belief propagation provides remarkable error cor¬ 
recting performance with reasonable time complexity. In recent 
days, it is easy to find practical applications of LDPC codes 
in wireless/wired communication systems and storage systems. 
Not only a practical importance but also recent theoretical 
advancement produces renewed interests in this field. Kudekar 
et. al E| proposed a new class of LDPC codes, that is called 
spatially coupled codes (SC-codes) and they provided theo¬ 
retical arguments on threshold saturation of SC-codes El- The 
origin of SC-codes is LDPC-convolutional codes that date back 
to the work due to Lelstrom and Zigangirov CD- Lentmaier et. 
al (6} showed an ensemble of an LDPC-convolutional code can 
have a higher threshold than that of a component LDPC code 
ensemble. Prom these works on SC-codes, it is unveiled that 
well-designed SC-codes have capacity achieving performance 
over symmetric memoryless channels. 

A burst erasure means a consecutive erased symbols. In 
many practical situations, we can observe occurrences of burst 
erasures due to slow fading in mobile wireless communication, 
buffer overflow at a congested router in a packet based 
network, and media flaw in a magnetic recording system. 
A strong erasure correcting code should have high erasure 
correcting capability not only for memoryless random erasures 
but also for burst erasures. Ohashi et. al 0 pointed out that 
SC-codes are not immune to burst erasures compared with 
conventional LDPC codes such as regular LDPC codes. In 
a typical decoding process of SC-codes, reliabilities of bit 
estimation gradually improves from both side into inside as 


a domino toppling. Since a burst erasure interferes the propa¬ 
gation of a wave of such reliable estimations, it causes severe 
degradation on decoding performance. In order to overcome 
this difficulties, they proposed a new class of multidimensional 
SC-codes that shows higher immunity against burst erasures. 

It is known that the burst erasure correcting capability of 
LDPC codes depends on a column order of parity check matri¬ 
ces of LDPC codes m- This is because the minimum length 
of stopping sets determining the burst correcting capability 
depends on the column order of a parity check matrix. In order 
to enhance the burst erasure correctability, several heuristic 
algorithms to improve the column order have been prensend 
by Wadayama 0, Paolini and Chiani E3, Hosoya et al. 0. 
Of course, the column order of a parity check matrix does 
not affect the decoding performance over memoryless erasure 
channels. 

In this paper, we will propose a class of column permu¬ 
tations that is called band splitting permutations suitable for 
( l,r,L ) SC-codes. A band splitting permutation is applied to 
the base matrix of ( l,r,L ) SC-codes having a single diagonal 
band and it results in a column-permuted base matrix with 
several diagonal bands. By lifting up the permuted base matrix, 
we can obtain a parity check matrix of a permuted (l, r, L) SC- 
codes. It will be proved that an appropriate band splitting per¬ 
mutation produces permuted (l, r, L) SC-codes that have near 
optimal minimum length of stopping sets. The permuted SC- 
codes constructed in such a way have burst erasure correcting 
superior to those of conventional SC-codes. Upper and lower 
bounds on the minimum length of stopping sets to be proved 
in this paper can provide theoretical performance guarantees 
for burst erasure correcting capability of permuted SC-codes. 

The outline of this paper is as follows. Section 2 provides 
notion and fundamental definitions required throughout this 
paper. Section 3 presents several theorems regarding stopping 
sets in a base matrix. The band splitting permutations will 
be defined and analyzed in Section 4. Results on computer 
experiments will be shown in Section 5. 

II. Preliminaries 
A. (I, r, L) SC-codes 

In this subsection, the definition (l,r,L) SC-codes pro¬ 
posed by Kudekar et al. 0 is reviewed. The (l,r,L) SC- 
codes belong to the class of protograph LDPC codes and 
its parity check matrix can be obtained by lifting up the 
base matrix B(l,r,L). The base matrix B(l,r,L) is a binary 
(L + l — 1) x kL matrix (k = r/l) and its structure is illustrated 


k 
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Fig. 1. Base matrix B(l,r,L) that defines ( l,r,L) SC-codes (gray area is 
filled with symbol one). 

in Fig. U] The parameters l and r represents the column weight 
and maximal row weight of B(l, r, L), respectively. We assume 
that the ratio k = r/l is integer throughout the paper. The 
parameter L denotes the number of sections. 

A parity check matrix of an (l,r,L) SC-code can be 
obtained by lifting up the base matrix B(l,r,L). A lift- 
up process is summarized as follows: For each element one 
in B(l,r,L), we can replace it with any binary M x M 
permutation matrix. The zeros in B(l, r, L) should be replaced 
with a binary M x M zero matrix. Let a parity check matrix 
obtained by the above process be H. The binary linear code 
defined by H is called an ( l,r,L ) SC-code. The size of the 
permutation matrices, M, is said to be the lift up factor. The 
number of rows of H is M(L + l — 1) and the number 
of columns is MkL. The design rate of ( l,r,L ) SC-codes, 
R(l,r,L), is thus given by 

R{1, r, L) = 1 — -j- — (1) 

B. Maximal correctable burst length 

Yang and Ryan fl2ll introduced a measure for burst erasure 
correcting capability of LDPC codes that is called the maximal 
correctable burst length. Let if be a parity check matrix 
that defines an LDPC code. The maximal correctable burst 
length of this code is denoted by W max (H). The meaning of 
W max (H) is the following. A burst erasure is a sequence of 
consecutive erasures occurred on an erasure channel. In this 
paper, we assume that only single burst erasure occurs in a 
code block. If the length of a single burst erasure is less than 
or equal to W morc (fT), it can be perfectly corrected by belief 
propagation (BP) decoding for erasure channels. On the other 
hand, there exists a single burst erasure of length W max (7/) + 1 
that cannot be corrected with erasure-BP. Namely, W ma x{H) 
represents the maximum guaranteed correctable length for any 
single burst erasure. As a related measure for burst erasure cor¬ 
recting capability, we here introduce the maximal correctable 
burst ratio defined by X max = W max (H)/n, where n is the 
code length. This quantity is useful for studying asymptotic 
behavior of the burst correcting capability. 

C. Stopping sets and maximal correctable burst length 

1) Stooping sets: Let H = (hi, b. 2 , ■ • •, h n ) € F™ xra be a 
parity check matrix. The vector hi is the i-th column vector of 
H. A sub-matrix of H consists of a subset of column vectors 
in H ; namely a sub-matrix of H has the form: 

H{ il _ iu} = (h il ,h i2 ,...,h iu )G F™ x “. (2) 


The subscript in Hu^ ...,i u } represents the column indices of 
H corresponding to the column vectors in the sub-matrix. 

Definition 1 (Stopping sets |fl3j ) Let II be a parity check 
matrix and S = {ii, 12 , ■. ■, i u } C [1, n] be an index set. The 
notation [a, b] denotes the set of consecutive integers from a 
to b. If the sub-matrix Hs has no rows with weight one, the 
index set S is said to be a stopping set. 

It is well known that stopping sets are closely related 
to correctability of erasure patterns if we exploit erasure-BP. 
Assume that a transmitted word is a codeword of the code 
defined by H and that some symbol erasures happen over the 
channel. Let E = {ei, e 2 ,..., e w j C [l,n] be the indices 
corresponding to the symbol erasures. This erasure pattern 
cannot be corrected with erasure-BP if there exists a non-empty 
stopping set S satisfying S C E. This fact indicates that the 
set of stopping set in H determine W max (H) lfl4l 

Assume that H = (hi, b- 2 ,..., h n ) G F™ xn is given and 
an index set S = {ii, ■ ■ •, i u } C [1, n] is given as well. The 
length of S, that is denoted by Len(S), is defined by 

Len(S) = 1 + max |a — b\. (3) 

a,b£S 

Let us denote the set of non-empty stopping sets of H by 
Q(H ). The span of H, Span(H), is defined by 

Span(H) = min Len(S). (4) 

ScQ(H) 

It is clear that a burst erasure of length shorter than Span(H) 
cannot cover any non-empty stopping set in H. This means 
that we have 

W max (H) = Span(H) - 1. (5) 

Note that the quantity W max (H) can be evaluated efficiently 
by using erasure-BP [fl2l . From the definition, we can see that 
Span(H) strongly depends on the order of the column vectors 
in II. It has been shown that an appropriate rearrangement of 
column order can increase the span of LDPC codes ®, a 

2 ) Irreducible stopping sets: We provide the definition of 
irreducible stopping sets that will be required in the next 
section. 

Definition 2 (Irreducible stopping sets) Let S C [l,n] be 

a non-empty stopping set of H. If removing any subset of 
elements from S yields an index set that is not a stopping set, 
then S is said to be irreducible stopping set. 


From the above definition of irreducible stopping sets, it is 
straightforward to see that the inequality Len(S') < Len(S) 
holds for a pair of nested stopping sets where S' is a stopping 
set and S' C S is an irreducible stopping set in S. From this 
inequality, we have 

Span(H) = min S ' G Q’(H)Len(S'), (6) 

where Q'(H) is the set of irreducible stopping set of H. This 
means that we only need to focus on the set of irreducible 
stopping sets when we discuss the span of H. 









III. Irreducible stoping sets in base matrix 

In this section, we will prepare several theorems regarding 
the maximal correctable burst length that are required for the 
argument in Section IIVI 

A. Maximal correctable burst length of base matrices 

Sridharan et. al ED studied the maximal correctable burst 
length of protograph LDPC codes. They showed a tight re¬ 
lationship between W max (B) and W max (H) where IT is a 
parity check matrix obtained by lifting up a base matrix B. 
The next theorem states this relationship. 

Theorem 1 (Maximal correctable burst length ( lH6l )) 

Assume that a base matrix B £ F™ xra is given. Let H be a 
parity check matrix obtained by lifting up B. The following 
inequalities hold: 

C W max {B) - 1 )M < W max (H) < ( W max (B) + 1 )M. (7) 

Theorem [T| indicates that the maximal correctable burst 
length of a protograph LDPC code is nearly determined by 
W max (B). This means that an appropriate column permutation 
for a base matrix B might be able to improve the maximal 
correctable burst length of a resulting photograph code. Of 
course, ( l,r,L) SC-codes belong to the class of protograph 
LDPC codes. It is reasonable to devise an appropriate column 
permutation for B(l, r, L), which will be discussed in the next 
section. 


B. Irreducible stopping sets in B(l,r,L) 

The maximal correctable burst length of the base matrix 
B(l,r, L) is determined by the set of irreducible stopping sets 
in B(l,r,L). In this subsection, we will show a structural 
property on the set of irreducible stopping sets in B(l,r,L). 

Let us denote the base matrix of the (7, r, L) SC-codes as 

B(l,r,L) = {b 1 ,b 2 ,...,b kL )£ ®T xfei . 

A block Ti(i £ [1, L]) that is a subset of indices is defined by 

Ti = {(i - 1 )k + 1, (i - 1 )k + 2, l)fc + k}. (8) 

From the structure of B(l,r,L) ( i.e ., Fig. [T). it is easy to 
see that b a = bp holds if and only if a. B £ Ti. The next 
theorem characterizes the structure of irreducible stopping sets 
in B(l , r, L). 

Theorem 2 (Irreducible stopping sets of base matrix) 

The set of irreducible stopping sets in the base matrix 
B(l,r, L) is given by 


0). Let us focus on the first nonzero element of the first column 
of the sub-matrix B(l,r, Due to the assumption 

that b a / bp and the definition of 11(1, r. L), it is evident 
that the row corresponding to the first nonzero element has 
Hamming weight 1. This means that S cannot be a stopping 
set in this case. 

By using this sufficient condition, we can immediately 
show that any stopping set of B(l, r, L) contains two different 
indices which belong to the same block. In other words, any 
stopping set must contain (a,/3) satisfying b a = bp(a (3). 
If a stopping set without such a pair exists, it contradicts the 
sufficient condition shown above. 

It is clear that a pair of indices (a, j8)(a, [3 £ [1, kL\, a ^ 
0) is an irreducible stopping set if both indices a and /3 belong 
to the same block. The last job is to show that there are no 
irreducible stopping sets with size larger than 2. Suppose that S 
is an irreducible stopping set with size larger than 2. From the 
above argument, S must contain at least a pair of two elements 
that belong to the same block. Since such a pair constitutes an 
irreducible stopping set, it contradicts the assumption that S is 
an irreducible stopping set. This completes the characterization 
of the set of irreducible stopping sets of B(l,r,L). □ 

C. Burst erasure correcting capability of ( l,r,L) SC-codes 

An immediate application of Theorem [2] is to analyze the 
burst erasure correcting capability of ( l,r,L ) SC-codes. The 
size of irreducible stopping set is two and the minimal length 
of the stopping set is thus two; we have Span(B(l , r, L)) = 2. 
This gives W max (B(l, r, L )) = 1 and we can utilize Theorem 
U] to obtain lower and upper bounds on maximal correctable 
burst length of ( l,r,L ) SC-codes: 


0 < W max (H) < 2 M, 


( 10 ) 


where H represents a parity check matrix of (l, r, L) SC-codes. 
By dividing both sides in (ITOt by the code length kLM, we 
have inequalities for the maximal correctable burst ratio: 


0 < 


< 


kL 


( 11 ) 


Q'(B(l, r, L)) = {{a, (3} \ a,/3 £ T t ,i £ [1, L]}. 


(9) 


It is clear that A max converges to zero when L goes to infinity. 
This inequality presents that the conventional ( l,r,L ) SC- 
codes have poor burst erasure correcting capability in the 
asymptotic regime when L —>• oo. This result justifies the 
observation made by Ohashi et. al (7). 

IV. Band Splitting Permutations 

In this section, we will propose band splitting permutations 
(BSP) for the base matrix B{l,r,L). The BSP is designed to 
improve the span of B(l,r,L). 


The theorem states that an irreducible stopping set consists of 
two column indices belonging to the same block. 

(Proof) Suppose an ordered index set S = {ji, ■ ■ ■ ,j u ) Q 
[l,n] is given where j\ < J 2 < ••• < ju■ The sub¬ 
matrix corresponding to S is written as B(l, r, = 

( bj x , bj 2 ,, bj u ). 

We will first show a sufficient condition that S is not a 
stopping set. Assume that b a bp holds for any a,/3 £ S(a ^ 


A. Definition 

When a BSP <Jk,L is applied to a base matrix B{1 , r, L ), we 
have permuted base matrix with multiple bands as shown in 
Fig. [2] The formal definition of BSP <Jk,L is given as follows: 
According to Cauchy’s two-line notation on a permutation, the 
permutation ak,L is described as 


Ofc.L = 


1 

/(l) 


2 

/( 2 ) 


kL 

f(kL ) 


( 12 ) 



Fig. 2. The structure of conventional base matrix and permuted base matrix. 


The second row of two line notation, i.e., the bijective function 
/ on [1, kL], is defined by 

(/(!) /(2) ••• f(kL)) = (a 1 a 2 ••• a k ), (13) 

where a,;(z £ [1, fc]) is given by 

ai = ( 1 1 + k 1 + 2 k 1 + (L —l)fc ) 


a k = ( k k+k k+2k fc+(L —l)fc ). (14) 


The permutation a k ,L can be seen as a block interleaver 
of interleaving depth k. Applying a ky L to the base matrix 
B(l, r, L) = bfci), we obtain a column permuted 

version of a base matrix B*(l,r,L ) = , &/(*,£)). 

For example, when k = 2, L = 3, we have 


<+ = ( 1 3 5 ), a 2 = 

( 2 4 

6 ) 


(15) 

( 

1 2 

3 

4 5 6 

\ 



ct 2,3 = 1 

1 3 

5 

2 4 6 



(16) 

Applying ct 2 ,3 to B{ 3,6,3), 

the 

permuted base 

matrix is 

obtained as 








/ 1 

0 

0 1 0 

0 ^ 




1 

1 

0 1 1 

0 



£*(3,6,3) = 

1 

1 

1 1 1 

1 


(17) 


0 

1 

1 0 1 

1 




^ 0 

0 

1 0 0 

1 ) 




Let us define k x L matrix A by 
/ a x 

A= : 

\ afc 


( 18 ) 


It is easy to see that z-th column of A corresponds the block 
Ti which is defined by ©. This implies that column vectors 
in B{l,r,L) belonging to the same block are rearranged in 
B*(l,r,L) as apart as possible. This property enhances the 
span of the base matrix. 


B. Bounds on maximal correctable burst length 

By lifting up the permuted base matrix B*(l , r, L ), we can 
obtain a parity check matrix of permuted (l,r,L) SC-codes. 
The next theorem provides upper and lower bounds on the 
maximal correctable burst length of permuted ( l,r,L ) SC- 
codes. This is the main contribution of this work. 


Let IT be a parity check matrix obtained by lifting up 
B*(l , r, L) with the lift up factor M. The maximal correctable 
burst length W max {H) of the permuted SC-code satisfies the 
following inequalities: 

(L — l)M < W max (H) < (L + l)M. (19) 

(Proof) Assume that S = {ii, i 2 , ■ ■ ■, i u } Q [1, rz] is a 
stopping set of B{l,r, L). The BSP maps S to 

Note that S* is also a stopping set of B*(l,r,L ) because 
B*(l, r, -^){/-i(i 1 ),/-i(i 2 ),...,/-!(»„)> contains a row of weight 
1 as well. This means that there is one-to-one correspondence 
between stopping sets in B{l,r,L) and those in B*(l,r,L). 
Theorem U indicates that a non-empty irreducible stopping set 
consists of two indices in the same block. Assume that a pair 
a, f3 £ [1 ,kL\ is such a pair of indices. From a definition 
of the matrix A in (IT8l) . it is clear that a and /3 belong 
to the same column in A. The definition of / in ( fTTl i thus 
leads to the inequality |/ -1 (a) — / -1 (/3)| > L that implies 
the length of irreducible stopping sets in B*(l,r,L) is larger 
than or equal to L + 1. Note that the equality holds when a 
and /3 are consecutive. From the definition of the span ©, 
we thus have Span(B*(l,r, L)) = L + 1 and this implies 
W m ax(B*(l,r, L)) = L. By using Theorem Q] the claim of 
this theorem is obtained. □ 

The inequalities of Theorem 0 indicates that the maximal 
correctable burst length of the permuted (l,r,L) SC-codes is 
proportional to the number of sections L. The inequality ( 1 1 () b 
indicates that the maximal correctable burst length does not 
depend on L for the case of the conventional (7, r. L) SC- 
codes. This result clearly shows the advantage of the permuted 
SC-codes over the conventional {i.e., non-permuted) SC-codes 
with respect to the maximal correctable burst length. 


C. Maximal correctable burst ratio 


In this subsection, we focus on the maximal correctable 
burst ratio X max of the permuted (l,r,L) SC-codes. 


By dividing both sides in ( IT9l > by the code length kLM, we 
can obtain following inequalities for the maximal correctable 
burst ratio: 


L - 1 
kL 


+ A 


< 


kL 


( 20 ) 


From (l20b . it is clear that \ max converges to 1/fc when L —>• 
oo. On the other hand, the design rate f?(Z, r, L) of the (Z, r, L) 
SC-codes converges to 1 — 1/fc as L goes to infinity. From these 
results, we have 


lim {Xmax + R{l,r,L)) = 1 (21) 

L—too 

that indicates asymptotic optimality of permuted ( l,r,L ) SC- 
codes in terms of burst erasure correction with erasure-BP. 
Note that no binary linear code of length n with design rate r 
can correct burst erasures of length larger than n( 1 — r). 


V. Numerical results 


Theorem 3 (Bounds on maximal correctable burst length) We have seen that the maximal correctable burst ratio of 

Let B*{l,r,L ) be the permuted base matrix defined above. permuted SC-codes can be approximated by X rnax ~ 1/fc 
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Fig. 3. Relation between L and A max of (l,r,L) = (3,6 , L) SC-codes 


Fig. 4. Histogram of A max of ( l,r,L ) = (3,6,32) SC-codes (M = 40, 

1000 samples) 


when L is large enough. We will here show the relationship 
between A max and L when L is finite. Figure 0 presents the 
bounds on A max and the BP threshold of (3, 6,L) SC-codes. 
The horizontal axis represents the number of sections L and 
the vertical axis is the maximal correctable burst ratio A max . 
It can be observed that A max of the conventional (3, 6, L ) SC- 
codes is decreasing as L increases. On the other hand, we can 
see that A max of the permuted (3,6 ,L) SC-codes converges 
to 1/2. It should be remarked that A max of the permuted SC- 
codes is higher than the BP threshold 0(3, 6, L) when L > 80. 
For example, when L = 128, A max of the permuted SC-codes 
is 0.496 but the BP threshold 0(3,6,128) is 0.488. Assume the 
case where the lift up factor M —>• oo. A combination of the 
conventional (l, r, L) SC-codes and an ideal symbol interleaver 
that can convert a single burst erasure into a memoryless 
random erasures may achieve the A max = 0(3,6, L). Thus, 
the permuted SC-codes yields better asymptotic burst erasure 
correcting performance when L is large enough. 

Figure [4] presents histograms of X rnax of randomly per¬ 
muted (3,6,32) SC-codes and permuted (3,6,32) SC-codes 
( i.e ., proposed codes), which are obtained by computer ex¬ 
periments with 1000 samples for each. The lift up factor is 
assumed to be M = 40. A randomly permuted SC-code is 
generated as follows: a parity check matrix of conventional 
SC-codes is at first produced and a uniformly random column 
permutation is then applied to it. It can be observed that A mra ., : ’s 
of randomly permuted SC-codes are far less than those of the 
proposed SC-codes. This result suggests that it is not trivial to 
find a superior permutation which provides better burst erasure 
correcting capability than that of a systematically designed 
BSP. 
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